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A  NUMERICAL  MODEL  SIMULATION  OF  A  SUMMER  REVERSAL  OF  THE  BEAUFORT  GYRE 

Ruth  H.  Frailer 
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Paaela  G.  Posey 
Berkeley  Research  Associates 


Abstract.  Ice  drift  derived  froa  the  Polar 
Ice  Prediction  Systea  (PIPS)  is  used  to  study 
the  sea  ice  circulation  of  the  Beaufort  and 
Chukchi  seas  in  the  late  suaaer.  Model 
slaulatlons  for  the  years  1983,  1986  and  1987 
shov  that  the  anticyclonic  Beaufort  Gyre 
reverses  to  cyclonic  circulation  in  the  late 
suaaer.  This  result  is  due  to  a  change  in  the 
aean  ataospheric  surface  pressure.  Drifting 
buoys  froa  1983  confira  this  reversal  of  the 
Beaufort  Gyre. 

Introduction 

The  Polar  Ice  Prediction  Systea  (PIPS)  is  a 
sea  ice  forecasting  systea  designed  for  the 
Arctic,  vhich  has  the  capability  of  predicting 
ice  drift,  ice  thickness  and  ice  concentration. 
The  PIPS  aodel,  developed  jointly  by  the 
Naval  Ocean  Research  and  Developaent  Activity 
(NORDA)  and  the  Fleet  Numerical  Oceanography 
Center  (FNOC),  is  run  daily  by  FNOC. 

The  Hlbler  ice  aodel  (Hibler,  1979]  serves  as 
the  basis  for  PIPS.  The  ice  aodel  is  driven  by 
|atmospheric  forcing  froa  the  Naval  Operational 
'Global  Ataospheric  Prediction  Systea  (NOGAPS) 
(Preller,  1985]  and  aonthly  aean  ocean  currents 
and  ocean  heat  fluxes  derived  froa  the  Hibler 
and  Bryan  ice-ocean  aodel  [Hibler  and  Bryan, 
1987].  The  systea  forecasts  over  the  entire 
Arctic  basin,  the  Barents  sea  and  the  northern 
half  of  the  Greenland  and  Norveglan  seas. 

PIPS  is  run  every  day  aaklng  a  120  hour 
forecast  of  ice  drift,  ice  thickness  and  ice 
concentration.  Each  day  the  aodel  is  initial¬ 
ised  froa  the  previous  day's  24  hour  forecast. 
Once  per  veek,  the  aodel  derived  ice 
concentration  field  is  updated  by  Arctic  ice 
concentration  data  made  available  by  the  Naval 
Polar  Oceanography  Center  (NPOC).  No  additional 
observational  data  is  incorporated  into  the 
model  at  any  other  tine.  Along  vith  the  120 
hour  forecast  aade  by  FNOC,  a  24  hour  forecast 
has  been  made  and  archived  dally  by  NORDA,  since 
June  1986,  as  a  means  of  monitoring  the  model. 

This  study  vlll  focus  on  ice  drift  fields 
vhich  have  been  derived  by  PIPS  during  1987  and 
vith  ice  drift  derived  froa  the  Hibler  aodel 
tested  vith  1983  and  1986  NOGAPS  data.  The  only 
difference  betveen  the  PIPS  results  and  the  test 
cases  for  1983  and  1986,  is  that  the  PIPS  ice 
concentration  field  vas  updated  veekly  and 
resulted  in  a  better  estimate  of  the  ice  edge 
location.  This  difference  does  not  impact  the 
results  presented  in  this  study. 

The  ice  model  has  been  shorn  by  Hibler  [1979] 
to  be  capable  of  reproducing  the  average  annual 
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ice  drift  in  the  Arctic.  This  average  pattern 
of  ice  drift  is  defined  by  a  large  anticyclonic 
circulation  in  the  vicinity  of  the  Beaufort  sea 
called  the  Beaufort  Gyre  and  a  transpolar  drift 
vhich  begins  near  the  Siberian  coast  and  extends 
through  the  Fran  Strait  (Gordienko,  1958]. 

Valsh  et  al.  [1985]  found  that  the  Hibler  ice 
model,  vhen  driven  by  30  years  of  atmospheric 
forcing,  shoved  interannual  variability  in  the 
annual  mean  pattern.  The  Beaufort  Gyre  and  the 
transpolar  drift  existed  in  all  cases,  but 
varied  in  intensity  from  year  to  year.  They 
also  found  that  the  Beaufort  Gyre  varied 
seasonally  in  size  and  jhape.  Model  results 
averaged  over  30  years  shoved  that  the  gyre  vas 
large  during  spring  and  fall,  covering  the 
vestern  third  of  the  Arctic  basin.  In  vlnter, 
the  gyre  covered  the  vestern  quarter  of  the 
Arctic  and  in  summer,  it  vas  reduced  to  a  small 
gyre  in  the  Beaufort  sea. 

This  study  vill  shov  that  vhen  model  results 
are  further  broken  dovn  into  aonthly  means 
during  the  summer,  a  distinct  change  takes 
place  in  the  standard  pattern  of  ice  drift.  The 
Beaufort  Gyre  is  seen  to  reverse  from  anti¬ 
cyclonic  to  cyclonic  circulation  in  the  late 
summer.  This  reversal  has  also  been  observed  in 
Arctic  buoy  data  [McLaren  et.  al,  1987], 

The  Model 

The  Polar  Ice  Prediction  System  is  based  on 
the  dynamic/ thermodynamic  Hibler  ice  model 
[Hibler,  1979].  The  PIPS  ice  model  has  five 
major  components:  a  momentum  balance,  a 
viscous-plastic  ice  rheology,  an  ice  thickness 
distribution,  an  ice  strength  relationship  and 
an  alr-lce-ocean  heat  balance  [Preller,  1985]. 
This  study  is  mainly  concerned  vith  the  momentum 
balance  used  to  determine  ice  drift: 

mDu-mfkxu+T.  +T„-mg  grad  H  +  &  (1) 

Dt  a  v 

vhere  m  is  the  ice  mass  per  unit  area,  u  is^the 
ice  velocity,  f  is  the  Coriolis  parameter,  t 
and  t  are  the  air  and  vater  stresses,  g  is  The 
acceleration  of  gravity,  H  is  the  sea  surface 
dynamic  height  and  F  is  the  force  due  to  the 
variation  in  internal  ice  stress.  Ice  moves  in 
a  tvo  dimensional  field  vith  forcing  applied 
through  simple  boundary  layer  formulations.  The 
air  and  vater  stresses  are  defined 

T  -  P.C.J&J  (tf  cos  6+  k  x  U„  sind)  (2) 

T  -  P  c  (j  -  U  |  [($  -U)  cos  e  +  k  X  (ft  -  U) 

V  W  W  W  Xin  6]  (3) 

vhere  J  is  the  ice  drift,  ti  is  the  geostrophic 

vind,  U  is  the  geostrophlc*ocean  current,  C 

Y  ® 


89  2  21  03 


Pr«ll«r  and  Posey i  Numerical  Siaulation  of  the  Beaufort  Gyre  Reversal 


70 


ICEVaoCfflES 

- 


JULY  ICE  VELOCITIES 


AUGUST  ICE  VELOCITIES 


SEPTEMBER 


Accession  For  j| 

NTIS 

DTIC 

Unann 

Justl 

GRAAI  rjf 

TAB  fl 
ounced  Q 
Stoat Ion 

Bv 

Distribution/ 

Aval 

1st 

M 

lability  Codes 

Aval*!  and/or 
Special 

f 

at 

20  m/stci  , 


*"•  j  /•;  %'r  'amftml 


Fig.  1.  1983  aonthly  aean  nodel  ice  drift  from  a)  July,  b)  August  and  c)  September. 
Vector  scale  is  30  ca/sec.  1983  aonthly  aean  NOGAPS  geos  trophic  winds  for  d)  July, 
e)  August  and  f)  Septeaber.  Vector  scale  is  20  n/sec. 


and  C  are  the  air  and  vater  drag  coefficients, 
p  anS  p  are  the  air  and  vater  densities  and  9 
aftd  8  are  the  air  and  vater  turning  angles. 
Surface  pressure  fields  froa  the  atmospheric 
aodel  are  used  to  calculate  geostrophic  vind. 
Surface  stresses  are  then  estimated  from  the 
geostrophic  vinds  using  a  turning  angle  of  23° 
and  a  drag  coefficient  of  .8  x  10  .  This  drag 

coefficient  vas  chosen  to  give  the  best  agree¬ 
ment  betveen  PIPS  ice  drift  magnitudes  and  the 
observed  ice  drift  magnitudes  from  Arctic  buoys. 
For  the  oceanic  boundary  layer,  the  ice  velocity 
relative  to  the  currents  beneath  the  boundary 
layer  is  used  to  calculate  a  quadratic  vater 
drag.  The  turning  angle  and  drag  coefficient  - 
used  in  the  calculations  are  23°  and  5.5  x  10”J 
respectively  (HcPhee,  1980J. 

The  PIPS  domain  covers  the  central  Arctic, 
Barents  sea  and  the  northern  Greenland  sea. 

The  nodel  grid  spacing  is  127  km  and  is  defined 
on  the  northern  hemisphere  polar-stereographic 
grid  used  by  FNOC.  The  model  uses  a  6  hour  time 
step  in  order  to  resolve  daily  variability  in 
the  atmospheric  heat  fluxes. 


Results 

The  Hibler  ice  aodel,  vhen  forced  by  N0GAPS 
atmospheric  forcing  for  1983,  1986  and  1987, 
along  vlth  the  ocean  currents  and  heat  fluxes 
from  the  Hibler  and  Bryan  ice-ocean  model,  shows 
annual  and  seasonal  mean  ice  drift  similar  to 
that  found  by  Hibler  [1979]  and  Valsh  et  al. 
(1985].  An  interesting  trend  is  seen  vhen  the 
summer  season  is  broken  down  into  monthly  means. 
Vhen  forced  by  1983  N0GAPS  data,  the  ice  model 
shovs  a  small  antlcyclonic  circulation  confined 
to  the  Beaufort  sea  during  the  months  of  July 
and  August  (Figures  la  and  lb).  However,  the 
September  monthly  mean  ice  drift  (Figure  lc) 
shovs  a  distinct  cyclonic  circulation  covering 
the  Beaufort  and  Chukchi  seas  and  extending  into 
the  central  Arctic.  The  gyre  veakena  froa  July 


through  August  and  then  completely  reverses  froa 
its  normal  antlcyclonic  to  cyclonic  drift  in 
September.  During  October,  the  monthly  mean 
circulation  reverts  back  into  the  small 
antlcyclonic  gyre  seen  during  July  and  August. 

The  circulation  of  the  gyre  closely 
corresponds  to  the  monthly  mean  vind  forcing, 
Figures  1  d-f.  Geostrophic  vinds  from  N0GAPS 
indicate  mean  clockwise  motion  in  July.  In 
August,  the  clockwise  circulation  has  become 
very  weak  in  the  Beaufort  sea.  By  September, 
the  vinds  are  rotating  counterclockwise  over  the 
Beaufort  and  Chukchi  seas. 

During  1983,  the  Arctic  Ocean  Buoy  Program, 
run  by  the  University  of  Washington's  Polar 
Science  Center,  had  from  6-13  buoys  in  the 
Arctic  during  the  July  -  September  time  period 
[Colony  and  Munoz,  1985].  These  buoys  provide 
information  on  position,  temperature  and 
atmospheric  surface  pressure.  Colony  and  Munoz 
have  also  created  charts  shoving  dally  surface 
pressure  fields  over  the  Arctic.  During  the 
period  of  July  through  Septeaber,  1983,  these 
surface  pressure  charts  shov  a  change  in  the 
mean  pressure  in  the  Beaufort  sea.  During  July, 
high  pressure  dominated  the  Beaufort  sea  region 
approximately  67X  of  the  month,  vhile  in  August 
it  drops  to  about  551.  During  September  the 
average  surface  pressure  pattern  changes  and  lov 
pressure  dominates  the  Beaufort  sea  approxi¬ 
mately  701  of  the  month.  This  data  vas  compared 
to  maps  of  surface  pressure  fields  from  the 
National  Meteorological  Center  (NMC).  The  tvo 
data  sets  agree  very  veil  during  these  summer 
months.  Monthly  mean  surface  pressures  from 
Colony  and  Munoz  [1985]  (Figures  2  a-c)  also 
compare  very  closely  to  the  N0GAPS  aonthly  means 
(Figures  2  d-f).  One  vould  expect  good  agree¬ 
ment  betveen  these  tvo  data  sets  since  buoy  data 
are  incorporated  into  the  NOGAPS  initialization 
scheme.  Figure  2f  has  included  the  displacement 
of  four  different  buoys  in  the  region  of  the 
Beaufort  gyre  during  the  month  of  September. 

All  of  the  buoys  indicate  counter-clockvise 
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Fig.  2.  1983  aonthly  aean  surface  pressure  fields  derived  from  buoy  data  (Colony  and 

Munoa,  1985)  for  a)  July,  b)  August  and  c)  September.  1983  NOGAPS  monthly  mean 
surface  pressure  fields  for  d)  July  e)  August  and  f)  September.  The  displacement  of 
4  buoys  during  the  month  of  September  is  included.  Contour  interval  is  2  mb. 
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motion  during  this  month  in  agreement  with  the 
model  results  (Figure  lc). 

For  the  annual  average,  PIPS  ice  drift 
appears  to  be  an  almost  equal  combination  of 
both  the  annual  average  ocean  currents  and  the 
annual  average  wind  field.  Hovever,  on  daily, 
veekly  and  monthly  time  scales,  the  vlnd  field 
dominates  the  ice  motion.  Thus,  during  the 
summer  months  vhen  high  pressure  dominates  the 
Beaufort  sea  region,  an  anticydonic  gyre 
results.  Hovever,  vhen  a  lov  pressure  field 
dominates,  as  in  September,  1983,  the  monthly 
aean  circulation  of  the  gyre  is  cyclonic. 

This  summertime  reversal  of  the  gyre  is  also 
found  vhen  the  ice  model  is  driven  by  1986 
NOGAPS  atmospheric  forcing.  Hovever,  in  1986, 
the  reversal  appears  earlier  in  the  summer 
(Figures  3  a-c).  During  July,  anticydonic 
drift  appears  as  the  monthly  mean  in  the 
Beaufort  sea.  A  cyclonic  gyre  dominates  the 
Beaufort  and  Chukchi  seas  in  August.  In 
September,  the  cyclonic  gyre  still  exists  but 
has  shifted  »o  the  northeast.  As  in  the  1983 
case,  the  drift  responds  directly  to  the 
circulation  Induced  by  the  atmospheric  surface 
pressure.  The  1986  NOGAPS  surface  pressure 
fields  show  good  comparison,  on  a  day  by  day 
basis,  with  the  NMC  surface  pressure  fields. 
Lov  pressure  systems  dominate  the  Beaufort  and 
Chukchi  seas  during  August  and  most  of 
September.  The  dominance  of  lov  pressure 
observed  in  this  region  is  not  due  to  a 
stationary  system  but  rather  a  series  of  lov 
pressure  systems  that  move  into  and  out  of  the 
area.  The  longest  period  of  tine  that  a  lov 


pressure  system  was  observed  to  have  remained  in 
the  Beaufort  and  Chukchi  seas  vas  seven  days. 
High  pressure  systems  appeared  interspersed 
between  these  lows  but  remained  in  the  region 
for  a  shorter  period  of  time. 

Results  from  the  PIPS  model  during  1987 
(Figures  3  d-f)  show  a  large  anticydonic 
Beaufort  Gyre  in  July  which  weakens  dramatically 
during  August.  In  September,  a  reversed 
Beaufort  Gyre  exists  hovever  it  is  shifted 
to  the  vest  of  the  cyclonic  Beaufort  Gyre  found 
in  the  summers  of  1983  and  1986.  As  in  the  tvo 
previous  cases,  the  monthly  aean  drift 
correlates  closely  to  the  monthly  mean  vind 
forcing. 

Conclusions 

This  study  shovs  that  numerical  ice  models  do 
simulate  a  reversal  of  the  Beaufort  Gyre  in  the 
late  summer.  The  PIPS  model  vhen  driven  by 
monthly  mean  ocean  currents  and  heat  fluxes  and 
NOGAPS  atmospheric  forcing  shovs  a  similar  trend 
in  the  monthly  mean  ice  drift  during  the  three 
test  years;  1983,  1986  and  1987.  The  trend  is 
defined  as  a  deterioration  of  the  normal  anti- 
cyclonic  Beaufort  Gyre  of  spring  and  early 
summer  into  a  complete  reversal  of  the  gyre's 
circulation  in  August  or  September.  This  change 
in  the  gyre's  circulation  is  due  to  a  change  in 
the  basic  pattern  of  atmospheric  surface 
pressure  in  this  region  of  the  Arctic.  The 
area  is  normally  dominated  by  high  pressure 
systems.  Hovever  in  the  late  summer  lov 
pressure  systems  appear  more  frequently  in  the 
Beaufort  and  Chukchi  seas. 
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Pig.  3.  1986  Monthly  aean  model  ice  drift  from  a)  July,  b)  August  and  c)  September. 

1987  monthly  mean  PIPS  ice  drift  from  a)  July,  b)  August  and  c)  September.  Vector 
scale  is  30  cm/ sec. 


Similar  results  have  recently  been  reported 
by  McLaren  et  al.  [1987].  Upon  examining  seven 
years  vorth  of  buoy  data  (1979-1985)  from  the 
Arctic  Buoy  Program,  McLaren  et  al.  conclude 
that  there  is  "a  reoccuring  tendency  for 
reversal  of  the  aean  clockvise  gyre  in  late 
summer,”  due  to  a  change  from  anticydonlc  to  a 
quasl-stationary  cyclonic  circulation  in  the 
atmosphere  in  the  vicinity  of  the  Beaufort  gyre. 
They  also  state  that  previous  testing  of  the 
Hibler  model  by  Hlbler  and  Tucker  [1979]  and 
Valsh  et  al.  [1985]  "failed  to  indicate  any 
cyclonic  pattern"  in  a  summer  siaulation  of  ice 
notion.  This  study  shovs  that  the  Hibler  ice 
model  does  indeed  predict  a  reversal  of  the  gyre 
in  the  late  summer.  This  reversal  varies  inter- 
annually  being  stronger  during  1983  and  1986 
then  during  1987.  This  trend  nay  not  have 
appeared  in  the  work  previously  published  by 
Hibler  and  Tucker  and  by  Valsh  et  al.  because 
they  discussed  annual  and  seasonal  means.  The 
reversal,  vhich  appears  for  only  one  to  tvo 
months  in  the  late  summer,  is  not  seen  in  the 
annual  mean  and  nay  not  appear  in  the  summer 
seasonal  average.  Also,  the  appearance  of  this 
reversal  during  the  three  test  years  does  not 
prove  that  it  occurs  every  year. 
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